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Abstract 

We consider here an extension of a previous work by Neves and 
Serva, still unpublished, which estimates the amount of interbreeding 
between anatomically modern Africans and Ncandcrtals necessary for 
explaining the experimental fact that 1 to 4% of the DNA in non- 
African living humans is of Neandertal origin. In that work we con- 
sidered that Africans and Neandertals had the same fitness (neutral 
hypothesis) and Neandertal extinction was thus an event of fortune. 
In this work we consider that Africans had larger fitnesses. We show 
results for four values for the fitness difference: 1%, 5%, 10% and 20% 
and compare them with the corresponding neutral results. Some tech- 
nical differences with respect to the neutral case appear. We conclude 
that even with 1% fitness difference Neandertals extinction comes up 
in too small a time, so the neutral model looks more suitable for ex- 
plaining the known data on occupation of some caves in Israel for a 
very long time, alternately by Africans and Neandertals. 



*This work is supported by FAPEMIG. 



1 Introduction 



For many years the question of the origins of the Homo sapiens seemed 
to have been decided in favor of the Out of Africa model. According 
to this model our species descends from a small group of Africans 
which appeared circa 100 to 200 thousand years ago and eventually 
spread all over the world. In its expansion this group would have met 
other human groups which already lived at other parts of the world. 
The most well-known of these groups is Neandertals, which inhabited 
Europe and West Asia until circa 30 thousand years ago and then 
disappeared. According to the Out of Africa model, other human 
groups, Neandertals in particular, belonged to different species and 
did not contribute to the H. sapiens gene pool. A competing the- 
ory is the Multiregional Hypothesis [22], according to which humans 
are the result of parallel evolution of some more or less separate re- 
gional groups, with occasional interbreeding among them making sure 
a unique species resulted. 

Support for the Out of Africa model came from paleonthological 
and archaeological arguments. But the support given by genetics was 
decisive. In the 1980's [21 [23] it was shown that all living humans are 
very similar to each other from the point of view of mitochondrial DNA 
(mtDNA). So similar indeed that mtDNA data are compatible with all 
of us being descendants of a single woman, which was then called the 
mitochondrial Eve. The time and place in which the mitochondrial Eve 
lived can be inferred from the distribution of differences among living 
humans in mtDNA and from geographical correlations. According to 
the supporters of the Out of Africa model, if interbreeeding with other 
human groups, e.g. Neandertals, had occurred, then we would have to 
see some rather different mtDNA types among living humans. Those 
different types have never been experimentally found. 

More recent genetic data of another kind seemed to confirm this 
view. By 1997 sequencing technology had advanced to the point of be- 
ing able to study mtDNA extracted from a few Neandertal fossils [TT| 
[To] and it resulted that indeed Neandertal mtDNA is rather different 
of the types existing among living humans. As a result of such seem- 
ingly decisive arguments, the Multiregional Hypothesis was almost 
completely discredited, although some statistical results favoring it 
could still be found [21]. 

Maurizio Serva[T7|, [THl [T9] and myself in collaboration with Carlos 
H. C. Moreira[13l [12l [H] observed in independent works that the ge- 
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netic arguments based on mtDNA utilized to support the Out of Africa 
model were not sufficient to conclude that the Multiregional Hypoth- 
esis was wrong. In order to have conclusive data we would have to 
resort to nuclear DNA. Very recent experimental data showed that we 
were right. Green et al[8] managed in 2010 to produce a first draft of 
the Neandertal nuclear genome. By comparing this genome with that 
of living humans from different geographic origins, they found out that 
Neandertals are genetically closer to present day non-Africans of ei- 
ther origin than to Africans. Not only this discovery provided a direct 
proof of interbreeding among Neandertals and our African ancestors, 
but it also allowed to estimate that living non- Africans carry in their 
cell nuclei 1 to 4 % DNA of Neandertal origin. Moreover, the fact that 
all non- Africans seem to be equally closer to Neandertals suggests that 
the interbreeding might have occurred in the Middle East between 45 
to 80 thousand years ago, when both groups occupied that region and 
before Africans had spread to the rest of the world. The strict Out 
of Africa model was definitely disproved. A second blow came still in 
2010, when Reich et aljlG] discovered that another human group pre- 
viously known only through its distinctive mtDNA did also contribute 
to the nuclear genetic pool of living Melanesians. 

The natural question posed by these discoveries is how much fre- 
quent the interbreeding should have been in order to explain the 1 
to 4 % proportion of Neandertal DNA in living non- Africans. Myself 
and Serva joined our forces to give a first answer to this question in a 
still unpublished work[15j. In that paper we devised a neutral model 
able to calculate the probability distribution for the amount of Ne- 
andertal nuclear DNA in an African population as a function of the 
interbreeding frequency of that population with a Neandertal popula- 
tion. By neutral we mean that in our model Africans and Neandertals 
are supposed to have the same fitness. 

In this work we intend to analyze the same model dropping the 
neutrality hypothesis. As it will be seen, the present analysis is some 
senses simpler than the previous [15], because we are justified in ne- 
glecting statistical fluctuations, but it introduces some technical ques- 
tions rather different from the ones in our previous work. In the next 
section we will present the model introduced by Neves and Serva [T3] 
and give an accurate definition also of the nonneutrality hypothesis. 
At the same time, we will also review our previous results [15] in order 
to compare them with the ones we will present here. 
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2 The interbreeding model 



We work on the hypothesis that Neandertals and Africans hved to- 
gether for some time in a restricted area, probably in the Middle 
East. A natural assumption then is that, due to limitation of natural 
resources, human population in that area remained of constant size, 
say equal to N individuals. We will suppose that such population is 
divided into two subpopulations we label as 1 (Africans) and 2 (Ne- 
andertals). Moreover, we will suppose for simplicity that generations 
are non-overlapping and we will count generations from past to future. 
Reproduction is sexual and diploid and we consider random mating 
within each of the subpopulations. We also suppose that the subpop- 
ulations have lived isolated from each other for a long time. The time 
when subpopulations meet and start sharing the same environment is 
labeled as generation g = 0. At that time total population consisted 
then of two groups, each of which composed by individuals of a pure 
race, which would live in the same restricted area for some time. 

The numbers Ni{g) and N2{g) of individuals at generation g in 
each of the two subpopulations are not supposed to be constant, al- 
though their sum Ni{g) + N2{g) = N is. In general, individuals will 
randomly choose mates belonging to their own subpopulation. Genetic 
mixing among subpopulations will be accomplished by assuming that 
at each generation a number a of random individuals from subpopu- 
lation 1 migrates to subpopulation 2 and vice-versa the same number 
of random individuals from subpopulation 2 migrates to subpopula- 
tion 1. In other words, a pairs per generation are exchanged between 
subpopulations. Migrants will participate in the random mating pro- 
cess of their host subpopulations and contribute with their genes for 
the next generation just like any other individual in that subpopu- 
lation. Their offspring, if any, is considered as normal members of 
the host subpopulation. The parameter a introduced above may be 
non-integer and also less than 1. In such cases we interpret it as 
the average number of pairs of exchanged individuals per generation. 
Our hypotheses concerning exchange of individuals suggest that, al- 
though subpopulations shared the same environment, there was still 
some kind of cultural barrier which prevented them from freely mat- 
ing. Despite that, we suppose that there was no biological barrier for 
interbreeding. 

We will suppose that the relative fitnesses of Neandertals and 
Africans are respectively 1 and 1 + h, i.e. Africans have in average 
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1 + h children per each Neandertal child. Parameter /i > is thus the 
fitness difference between Africans and Neandertals. To be clear, an 
individual which was born among the Neandertal subpopulation, but 
was exchanged and had offspring among African subpopulation is con- 
sidered as African for what concerns fitness. Vice- versa African born 
individuals living among Neandertals will have a Neandertal fitness. 
This hypothesis implies that fitness is considered as a cultural tract, 
not a genetic tract. Such an assumption is obviously disputable, but 
we believe it is not at all absurd. Moreover it simplifies the mathe- 
matical treatment given here and should not produce many errors for 
the small values of a and h to be used here. 

Let x{g) = Ni{g)/N be the fraction of Africans in the population 
at generation g. We will define the history of the subpopulation sizes 
as the set of values of x{g) at all generations. In case the size of 
the population were infinite, then the subpopulation history would be 
determined by 

^ {l + h)x{g) ^ {l + h)x{g) 

' {l + h)x{g) + l{l-x{g)) 1 + hx{g) 

= Vhixig)). (1) 

For finite we have statistical fiuctuations and it is reasonable to 
suppose that the history will be specified by supposing that Ni{g + 1) 
is a binomially distributed random variable such that its expectation 
turns out to be Nrih{x{g)). More concretely, 

PToh{N,{g + l)=n,) = ^ ^Vhix{g)r{l-Vh{x{g) f-^' ■ (2) 

Of course, N2ig + 1) = N - Ni{g + I). 

In the nonneutral case /i > 0, it can be easily seen that and 1 are 
the only fixed points of r]h{x), being unstable and 1 stable. Thus, if 
were infinite, x{g) ^ 1 as g ^ oo. For finite N statistical fluctuations 
ensure that x{g) will be equal to 1 at some flnite g with very large 
probability if A^ is large enough, even for small h. It can be shown that 
the mean number of generations for Neandertal extinction diverges 
only with logA^, being thus rather insensitive to A^ in a large range 
of values. Moreover, for a fixed A^, the distribution of the extinction 
times is relatively narrow. Such features are illustrated in Fig. [T] 

On the other hand, the possible set of histories in the neutral case 
/i = is much more variable. The stochastic process for Ni{g) and 
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Figure 1: The plots show histograms for the time to extinction empirically 
obtained by 1000 simulations of the stochastic processes for the population 
history x, all with x(0) = 0.5. The two plots in the upper row are for the 
nonneutral case with h = 0.1 and the two plots in the lower row are for 
the neutral case h = 0. The plots in the left column refer to population 
size N = 100 and the ones in the right column are for = 1000. Notice 
that whereas the horizontal axes in the plots for h = 0.1 represent number 
of generations, in the plots for h = they represent number of generations 
divided by N. 
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N2{g) turns out to be the neutral Wright-Fisher model for two alleles 
at a single locusfT]- In this case, not only the time for extinction is 
random, but also which of the two subpopulations becomes extinct. 
Africans will be the extant subpopulation with probability x(0) and 
the mean number of generations until extinction [7J is 

- 2A^[x(0) Inx(O) + (1 - a;(0)) ln(l - x(0))] . (3) 

The extinction time for a fixed x(0) rather than being relatively in- 
sensitive to is in fact proportional to N . See the illustration also in 
Fig.[l} 

The expectation oi x{g+\) — x{g) is r]h{x{g))—x{g) in both neutral 
and nonneutral cases. Whereas this quantity is exactly zero in the 
neutral case, it is positive in the nonneutral case. This explains the 
large qualitative difference between the two cases. 

Our goal is to quantify the amount of Neandertal genes in Africans 
at the time either of Neandertal extinction or Neandertal assimilation 
by Africans. The technical definition of assimilation will be explained 
ahead when it appears. The fraction of Neandertal genes in Africans 
should be compared with the experimental estimate [8] that living non- 
Africans have 1 to 4% Neandertal DNA. 

A precise definition of Neandertal and African genes follows from 
our hypothesis that subpopulations have been isolated from each other 
for a long time before g = 0- This allows us to suppose that in many 
loci the two subpopulations will have different and characteristic alle- 
les. Therefore, we can assume that there exists a large set of alleles 
which are exclusive of Africans and the same for Neandertals. We will 
refer to these alleles respectively as African and Neandertal. At any 
generation g > any individual will be characterized by his/her frac- 
tions of African and Neandertal alleles. We define then yi{t) as the 
mean fraction of African alleles in African subpopulation at genera- 
tion g and 2/2 (i) o,s the mean fraction of African alleles in Neandertal 
subpopulation at generation g. The mean here is due to the fact that 
individuals in African subpopulation in general have different allelic 
fractions, but yi{t) is calculated by summing allelic fractions of all 
individuals in the African subpopulation and dividing by Ni{t). Sim- 
ilarly for 2/2 (i)- Of course at the initial generation we have 

yi(0) = l and 2/2(0) = . (4) 

Similar quantities might have been defined for Neandertal alleles, but 
they are easily related to yi (t) and 2/2 (t) and thus unnecessary. 
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The equations relating the mean allehc fractions at generation g+1 
with the mean aUeUc fractions at generation g were obtained in our 
previous workjl5j. In obtaining them we will further assume a mean 
field hypothesis, i.e. that the a individuals of the African subpopu- 
lation migrating to the Neandertal subpopulation all have an allelic 
fraction equal to yi{t) and vice- versa for Neandertal migrants. The 
above mean field assumption is a strong one and it is not strictly true. 
Nonetheless, it is a very good approximation if a is much smaller 
than l/log2 N. In fact, 1/a is the number of generations between two 
consecutive exchanges of individuals. The typical number of gener- 
ations for genetic homogeneization in a population of N individuals 
with diploid reproduction and random mating, calculated by Derrida, 
Manrubia and Zanette[ll O |6] and by Chang[3], is log2 A^. Thus, the 
condition that a is much smaller than 1 / log2 N makes sure that sub- 
populations are both rather homogeneous at the exchange times. 

The allelic fraction yi{g + 1) will be equal to yi{g) plus the con- 
tribution of African alleles due to the immigrating individuals of Ne- 
andertal subpopulation 2 and minus the loss of African alleles due 
to emigration. We remind that these loss and gain terms are both 
proportional to a and inversely proportional to the number Nx{t) of 
individuals in the African subpopulation. Similar considerations apply 
to 2/2(5 + !)• The equations below are then exactly the same as in our 
previous work[15J, suitably rewritten in order to take into account the 
similarities and differences between the neutral and nonneutral cases: 

^2(5 + 1) = jv(i-x(g)) ^i(g) + (1- jv(i"x(g)) ) y^ia) 

It is understood that the above difference equations must be supplied 
with initial conditions Q. 

Another condition we must impose on the above equations is that 
at any generation g we have both 

2a < Nx{g) and 2a < N{1 - x{g)) . (6) 

These conditions ensure that at generation g both subpopulations have 
more than 2a individuals, which implies that there exist individuals 
enough to be exchanged and also that after the exchange, migrants will 
not be the majority of any subpopulation. When the above conditions 
fail to hold we may consider that one of the subpopulations, even if 
not yet extinct, is so close to extinction that it has been assimilated 
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by the other subpopulation. In this case, there is no more sense to 
consider subpopulations as separated. 

It is possible to see that the non-assimilation conditions ([6]) imply 
that 

x{g){l-x{g))>^, (7) 

which will be used below to guarantee some nice mathematical prop- 
erties of the solutions to ([sj |4]) . 

In the neutral case, due to ^ it is natural to define a rescaled 
time as t = g/N. As a consequence, ([s]) become a system of ordinary 
differential equations [15j in the limit — )• oo, 

= -^)iyiit)-y2{t)) 

The same rescaling does not work in the nonneutral case, because it 
would lead to Neandertal extinction in rescaled time equal to zero. 

Another difference is that in the neutral shown by ([s]) , the 

relevant parameter for quantifying the amount of interbreeding is a, 
i.e. the number of exchanged pairs of individuals per generation. In 
particular, larger values of A^ imply more generations until one of the 
subpopulations is extinct, but the final values of the allelic fractions 
yi and y2 are independent of A^, depending only on a and on the 
particular realization of the history x{t). On the other hand, in the 
nonneutral case, the relevant interbreeding parameter is 

/3=^, (9) 

the fraction of the entire population exchanged per generation. 

A qualitative picture of the solutions of ([5| [4]) can be obtained in 
a way analogous to what we did before[15]. Introducing the auxiliary 
functions zi{g) = yi{g) - y2{g) and 22(5) = Viio) + 2/2(5), Eqs. ^ 
become 

= - x(9){l-x(g)) ^1(5) 




By using ([T]) , the first of these equations shows that the difference 
between yi and 7/2 is always positive and decreases in time until either 
Neandertals are extinct or assimilated. Moreover, this difference will 
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become close to zero if /? is large enough, or if the number of gener- 
ations until Neandertal extinction or assimilation is large. If this is 
plugged again in ([5|, rewritten as 

yi{g + 1) - yi{g) = -^^1(5) ^^^^ 
y2{g + 1) - y2ig) = T^ziig) 

we see that yi must decrease and y2 must increase with time. We also 
see that the rates of decrease of yi and increase of 2/2 are in general 
different. The transfer of Neandertal alleles to the African subpopu- 
lation and of African alleles to the Neandertals will be symmetrical 
only if x{g) = 1/2. Moreover, the transfer of African alleles to the 
Neandertal subpopulation and vice-versa is more effective at initial 
generations, when zi(g) is larger. 

The above qualitative view is enough for many purposes; for quan- 
titative purposes we may easily numerically iterate the difference equa- 
tions ([5]), thus obtaining yi and y2 at all generations. This is compu- 
tationally much easier than running simulations of the entire stochas- 
tic processes of reproduction and gene transfer, described thoroughly 
before [T5]. A good agreement between results of simulations of the en- 
tire stochastic processes and solutions to ([s]) has been illustrated [T5] 
and will not be repeated here. 



3 The deterministic approximation and 
results 

Let y( denote the final value of yi at the time at which the Neandertal 
subpopulation is either extinct or assimilated. In the neutral case|15j 
y( will depend on a and x(0) and, given x{0), will also depend very 
much on the particular realization of the history x{g) which follows. As 
we have no information on the value of x(0) and its subsequent history 
x{g) in the particular realization of the process which produced the 
H. sapiens population as it is nowadays, the best we could do in the 
neutral case was to produce many simulations with random values of 
x(0), random subsequent history and random a and look at y( in such 
simulations. After doing that we counted the number of realizations 
for a particular value of a such that y( belonged to the experimental 
interval[8j of 1 to 4%, i.e. 0.96 < y{ < 0.99. As a result of such a 
process we obtained[l5], in the neutral case, an empirical probability 
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density distribution for a such that y{ falls into the experimental 
interval. The result [15] is that the maximum of this distribution is 
for a ~ 0.013 and its mean value is a ~ 0.083. As the reader may 
notice, by taking random histories the only parameter which needs to 
be estimated in the neutral case is a. This small number of parameters 
to be adjusted was considered[lJ a strong point in our previous results. 

In the nonneutral case, unfortunately we have two extra parame- 
ters with respect to the neutral case: the fitness difference h and the 
total population size N. On the other hand we have an advantage: 
the fact that the expected value of x{g + 1) — x{g) is positive makes 
the dependence on the particular history x{g) much smaller than in 
the neutral case. As a consequence stochastic histories are very well 
approximated by the deterministic history, i.e. the history obtained 
in the N ^ oo limit, given by x{g + 1) = r]h{x{g)), supplemented by 
the stopping condition - without it, the history will never end - 

>min{l- 1,1-2/3} . (12) 

The above condition means that we stop the deterministic history 
when either there is less than one Neandertal individual, or Neander- 
tals were assimilated. We may then use the deterministic history in 
Eqs. ^ to calculate y(. Such a calculation will be referred to as the 
deterministic approximation for y(. 

In Fig. 1] we compare results from simulations of the complete 
stochastic processes of reproduction within subpopulations and mi- 
grations until Neandertal extinction or assimilation with the corre- 
sponding results in the deterministic approximation. The reader may 
see that agreement is very good. 

We may now proceed to our main objective, which is estimating 
the interbreeding parameter f3 for some chosen values of the fitness 
difference h: h = 0.2, h = 0.1, h = 0.05, h = 0.01. As the left part of 
Fig. [3] illustrates, and as remarked before, the dependence of y( on A^, 
all other parameters fixed, is very slight. As a consequence, we will 
produce results for a single value = 1000, which seems a reasonable 
value. For other values of N the results for /3 ~ but not for a - are 
approximately the same. 

The data in Fig. [3] also illustrate the range of values of allowed 
values of /3 for each x(0). By using the data in the figure we may 
calculate the empirical probability density distribution for the values 
of /3 such that y( lies in the experimental interval, a result analogous 
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Figure 2: In the left we show plots of a simulated history (dots) and the 
deterministic history (full line) for a population with size N = 200, initial 
fraction of Africans x(0) = 0.5 and fitness difference h = 0.2. In the right 
we show the histogram of the yl values obtained for 100 simulations of a 
population again with = 200, x(0) = 0.5 and h = 0.2. The value of 
the interbreeding parameter is /3 = 5 x 10~^, which corresponds to one pair 
of individuals exchanged at each 10 generations. For comparison sake, the 
mean of the presented y( data is 0.980636, while the corresponding value 
calculated in the deterministic approximation is 0.983291. 



to our main result [15J in the neutral case. This probabiUty density 
distribution is plotted at the right in Fig. [3] for h = 0.1 and = 1000. 
For this plot we have restricted the histories to those such that 

x(0) < 0.8 (13) 

as in our discussion[15j of the neutral case. The reason for this restric- 
tion is that larger values of x(0), as can be seen in the left part of Fig. 
[3| strongly shift the probability density for /3 to larger values. As the 
mean field hypothesis included in [5] is not accurate for large values of 
a, we add the assumption that Neandertals were at least 20% of the 
total population at g = 0. Apart from the mathematical reasons 



for the no Neandertal minority assumption (13), we remark that it is 
also reasonable from a historic point of view, since Neandertals and 
Africans seem to have coexisted for thousands of years in the Middle 
East. Of course the value 0.8 is somewhat arbitrary, but it is the 
same we have chosen [15] in the neutral case and is used here also for 
comparison sake. 

In the following figures we repeat the same calculations as in Fig. 
[3] to the values h = 0.05, h = 0.01 and finally h = 0.2. In all of them 



we use N = 1000 and the no Neandertal minority assumption (13). 
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Figure 3: At the left we show for h = 0.1 and values = 100, = 200 
and A^ = 1000 the regions in the plane x{0), (3 such that y{ lies in the 
experimental interval between 0.96 and 0.99. The picture was obtained by 
calculating y{ in the deterministic approximation for a grid of points in the 
plane and then drawing the level curves of y{ at levels 0.96 and 0.99. Notice 
that the regions for the different values of A^ are almost the same, mainly in 
the larger /3 region. At the right we plot the probability density distribution 
for the values of (5 such that y{ in the deterministic approximation lies in 
the experimental interval between 0.96 and 0.99, taking also into account the 
no Neandertal minority assumption (13). The parameters are h = 0.1 and 
A^ = 1000. The mean value of (3 is /3mean ~ 5.3 x 10"'^. The maximum is at 
/3^ax ~ 2.5 X 10-4 
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Figure 4: The same as in Fig. [3} but now h = 0.05 and we show only 
the region for = 1000. The mean value of /3 is /3mean ~ 2.7 x 10"''. The 
maximum is at /3max ~ 1-3 x 10~^ 
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xiO) 

Figure 6: The same as in Fig. |4| with h = 0.2 and = 1000. The mean 
value of (3 is /3mean ~ 9.3 X lO"''. The maximum is at /3max ~ 5.0 x 10~^. 

4 Discussion and conclusions 

We have presented here the extension to a nonneutral setting of a 
model [T^ originally intended to estimate the amount of interbreed- 
ing between anatomically modern Africans and Neandertals suppos- 
ing equal fitnesses for these groups. As commented along the text 
there are some subtle differences in the mathematical treatment of 
the two cases. Due to these differences the interbreeding parameter 
a introduced in our previous work was substituted hy /3 = a/N. In 
doing this, results turn out to be approximately independent of N, as 
illustrated in the left part of Fig. [Sj Comparison of the present results 
with the previous ones[T5| is possible only by choosing some value for 
A^. Although results in Figs. |3} |4j [5] and [6] were produced with the 
choice N = 1000, the A'-independence permits that we extrapolate 
them to other values for A^. 

In the left part of Fig. [7] we present a summary of the results 
of Figs. [3j [4j [5] and [6] along with the corresponding results for the 
neutral case h = 0. It turns out that, as expected, larger values 
for the fitness difference imply shorter times until Neandertals are 
extinct or assimilated and, as a result, interbreeding must be more 
intense in order that the necessary amount of Neandertals alleles be 
transferred to the African subpopulation. This explains why the values 
of f3 increase with h. 

Although we used /? instead of a in this work, we remind that the 
mean field hypothesis used in deriving Eqs. ([5| relies on the smallness 
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Figure 7: In the left picture we summarize the results of the neutral, h = 0, 
and nonneutral cases, giving as functions of h the positions /3max of the maxi- 
mum of the probability density distributions and /3mean of their means. Filled 
dots are data for /Smax and empty dots are data for /3mean- In the right picture 
we show for a total population = 1000 the number of generations necessary 
for extinction of Neandertals in the deterministic approximation as functions 
of the initial fraction of Africans x(0). The various curves correspond to the 
values h = 0.01, h = 0.05, h = 0.1 and h = 0.2. For fixed x(0) the number 
of generations is a decreasing function of h, which makes labelling of the 
curves unnecessary. For other values of N results do not change very much, 
as commented in the text. 
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of a. More exactly, the mean field hypothesis is accurate for a <C 
l/log2iV. For N = 1000, this means a <C 0.01, or /3 < 10"^. Many 
of the values found for /3 in Fig. [7] are larger than that. These values 
should be taken with care if A'^ = 1000, but if N were smaller, say 
N = 100, they may be considered accurate. We do not know any 
estimate for in the Middle East, but 10,000 individuals [20j has 
been suggested as a good estimate for the total human population in 
ancient times. 

Another question to be considered in applying the present results, 
not yet addressed, is the time necessary for the extinction of Nean- 
dertals. We show in the right part of Fig. [7] a plot of the number of 
generations until Neandertal extinction, in the deterministic approx- 
imation, as a function of the initial fraction x(0) of Africans for the 
various values for h used here. Taking one generation to be roughly 20 
years, we see that the times range from 600 to 1000 years for h = 0.2 
to 11,000 to 20,000 years for h = 0.01. 

According to Bar-Yosef[9j, occupation of the Middle East by Ne- 
andertals and Africans can be compared with a long football game, 
in which teams alternated in their ability to dominate the game field. 
Such a situation is more likely to be described by the neutral model, 
but a nonneutral situation, in which Africans had a slight advantage, 
could be possible. Our present results show that even a 1% fitness 
advantage in favor of Africans, the smallest h value we investigated, 
seems to be too large, as the caves of Skuhl and Kafzeh, in Israel, 
alternated between Africans and Neandertals several times over a pe- 
riod of more than 130,000 years. For smaller values of h it is possible 
that statistical fluctuations make the deterministic approximation em- 
ployed here useless. As our previous work on the neutral model[15j did 
take into account the unavoidable statistical fluctuations, we believe 
- up to now - that the neutral model is more suitable for explaining 
both the amount of Neandertal DNA in present day humans and the 
long time Africans and Neandertals lived in the Middle East. 
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